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O ■ Abstract 

O ' We present a generalization of the five dimensional multigravity models to six dimen- 

sions. The key characteristic of these constructions is that that we obtain solutions which 
do not have any negative tension branes while at the same time the branes are kept flat. 
A ■ This is due to the fact that in six dimensions the internal space is not trivial and its curva- 

ture allows bounce configurations with the above feature. These constructions give for the 

first time a theoretically and phenomenologically viable realization of multigravity. 
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1 Introduction 

Over the past two years there has been considerable activity studying brane configurations 
in which gravity is modified in large scales motivated by higher dimensional construc- 
tions |]-§|. The interesting feature of these scenarios is that modification of gravity at 
cosmological scales was triggered by microphysics considerations which have to do with the 
structure of the extra dimensional space of usually Planck scale volume. These scenar- 
ios are in principle testable by the CMB power spectrum JF| and gravitational lensing || 
observations. 

The prototype multigravity scenario was the "H V" bigravity model where, apart of 

the graviton zero mode, an anomalously light graviton Kaluza-Klein (KK) state is present. 
This gave the possibility that gravity at intermediate scales (1mm < r < 10 26 cm) is me- 
diated by both the massless graviton and the ultralight state, leading to the terminology 
"bigravity" . At large scales were the Yukawa suppression of the massive state turns on, 
the gravitational constant decreases to a value which for example in the symmetric config- 
uration is half of its value at intermediate scales. For asymmetric configurations one can 
have a situation were the the massive graviton dominates and gravity switches off at large 
scales or, on the other hand, where the massless mode dominates and gravity does not 
change at all. Similar constructions of freely moving negative tension branes can also give 
qualitatively different modifications of gravity at large scales. The quasi-localized gravity 
model of Gregory, Rubakov and Sibiryakov || has no normalizable zero mode but gravity 
at intermediate distances is generated by a "resonance" -like ||,[nj coupling of the low part 
of the continuum of the KK spectrum. In this construction, the gravitational law at large 
scales changes from four dimensional to five dimensional. The change of the nature of the 
gravitational law was observed also in the the crystalline universe model || where gravity 
at intermediate distances is generated by the first low lying band of KK states and at large 
scales gives again the five dimensional law of gravity. 

So far our discussion has considered the gravitational force due to tensorial gravita- 
tional modes. However, in the gravity sector there exists a moduli corresponding to the 
fluctuation(s) of the negative tension brane(s), the radion, which changes the large scale 



phenomenology of these models. In particular, the radion is a ghost field fTT|^I3|, i.e. it has 



wrong sign kinetic term, so at large scales it will dominate in all the above constructions 



giving scalar antigravity Jl4j instead of a reduction of the Newton's constant or a change 



of the nature of the gravitational law. However, the very fact that there is a physical ghost 



in the spectrum makes these models unacceptable. Note that the essential characteristic of 
all these models is the bounce of the warp factor generated by the presence of the moving 
negative tension branes. This association between the bounce in the warp factor and pres- 
ence of negative tension branes is an unavoidable feature in flat (five dimensional) brane 
models and is linked to the fact that the weaker energy condition is violated at the position 
of the bounce |15, 16|. As a result one must have the unacceptable ghost field. 



A way out of this difficulty is to consider AdS^ branes instead of flat branes (see [17] for 
a different possibility involving an external four-form field). In this case the warp factor has 
a " cosh" form which naturally generates a bounce without any need of floating negative 
tension branes. Indeed, it is straightforward to replicate the bigravity and the crystalline 
model in this framework. The bigravity model is converted to a " + +" model |]TB[ and 
the crystal model to an infinite array of "+" branes. The quasi-localized model on the 
other hand cannot be reproduced because the presence of the cosmological constant on 
the branes prevents the warp factor from being asymptotically constant. By sending the 



second brane in the "++" model to infinity, one obtains the locally localized model of fl"9|.[20 
where gravity is mediated only by the ultralight state since the zero mode becomes non- 
normalizable. These models, even though they solve the theoretical difficulty of the ghost 
radion(s) in the flat brane models, face phenomenological difficulties. The presence of a 
remnant negative cosmological constant is at odds with observations, and furthermore it 
turns out that all large scale modifications of gravity predicted by these models, are hidden 
behind the AdS horizon. 

It is interesting to note at this point how these models confront the relativistic predic- 
tions of four dimensional Einsteinian gravity. It has been known for a long time that the 
structure of the massive graviton propagator is different form the massless one in flat space, 
something that is known as the van Dam - Veltman - Zakharov discontinuity |21|,|22[ . This is 



due to the extra polarization states of the massive graviton and in particular the scalar-like 
mode that does not decouple in the limit of vanishing mass. As a result, gravity generated 
by a massive graviton case predicts a 25% discrepancy in the bending of light observation 
if one insists in the validity of Newton's law for static bodies or vice versa. This apparently 
rules out all models of gravity based on a massive graviton [|I(J] since the observed bend- 



ing of light by the Sun agrees with the General Relativistic result to an accuracy of 1% 
(see also [^3|] for a discussion on the relation of this observation with physics arising from 
higher curvature terms in four dimensions). In the models of the flat branes, the ghost 
radion(s) at intermediate distances cancel the troublesome additional polarization state(s) 



of the massive graviton(s), restoring the phenomenologically favoured tensorial structure of 
the graviton propagator fl2|,|14l,|2"4||. In the models with AdS branes it was shown that the 
extra polarization state is practically decoupled due to an unusual feature of the graviton 
propagator structure in AdS. In more detail, it was shown p5|-p8| that the van Dam - 
Veltman - Zakharov no-go theorem can be evaded at tree level in (A)dS in the case where 
m/H — > 0, where m is the mass of the graviton and H the "Hubble" of the AdS space, 
a precise realization of which was the " + +" brane model. Thus, in the AdS branes no 
phenomenological constraint from the bending of light experiment is applicable provided 
this condition holds. The discontinuity in the pure four dimensional theory with a massive 



graviton reappears at the quantum level [29. 30] but one should keep in mind that purely 



four dimensional theory of massive gravity is not well defined J5I]. If instead one starts 
from the higher dimensional theory in which the mass terms appear dynamically and if one 
takes into account the whole tower of the KK states, this quantum discontinuity is expected 



to be absent, as pointed out in [13]. Let us also note that in a recent paper [|32| it was 



suggested that the extra polarizations of the massive state in the one brane model can be 
gauged away and thus the discontinuity is absent. 

More generally, the argument of van Dam - Veltman - Zakharov that no mass for the 
graviton is allowed in a flat background is not entirely robust. A long time ago Vainshtein 



33] showed that if one considers the bending of light, the lowest order approximation is 



not valid because the calculation does not take into account the characteristic mass scale of 



the problem. A more detailed non-perturbative analysis presented in the recent paper [|34 
supports the original Vainshtein argument that there is continuity in the gravitational 
potential as the mass of the graviton goes to zero for distances smaller that a critical 
one which, in the case of the Sun, is of order of the size of the solar system. The linear 
approximation is valid for distances greater than this critical one, where the discontinuity 
reappears. Thus the observation of light bending by the Sun does not rule out the massive 
graviton proposal. 

In the light of this encouraging result it seems appropriate to look for multigravity 
models of flat branes without ghost fields. As we have discussed, in five dimensions it is 
impossible to have flat brane multigravity models without negative tension branes. For this 
reason, we will consider models in six dimensions. The literature on six dimensional con- 



structions is already quite rich |35Tp3|. In this paper we explicitly show that the difficulties 
of the five dimensional models can be evaded in six dimensions. It is possible to construct 
multigravity models with only flat positive tension branes. The branes which localize grav- 



ity in this setup are four-branes, but one of their dimensions is compact, unwarped and 
of Planck length. Thus, in the low energy limit the spacetime on the brane appears three 
dimensional. In order that these constructions are realized it is crucial that the tensions 
and/or the bulk cosmological constant are inhomogeneous. The five dimensional construc- 
tions in this setup come into two types. One of them involves conical singularities at finite 
distance from the four branes. These conical singularities support three-branes^ which can 
be of positive tension if one has an angle deficit, of zero tension if one has no angle deficit 
and of negative tension if one has an angle excess. The other type of the multigravity 
models has no conical singularity at all. It is crucial to note that there is no five dimen- 
sional effective theory for these constructions, otherwise one would get all the problems 
faced in the five dimensional constructions. The low energy effective theory is directly four 
dimensional. 

The organization of the paper is as following. In section two we present the cylindrically 
symmetric single four-brane warped models and show the cases where gravity is localized on 
the four-branes. In section three we paste two of the single four-brane solutions and obtain 
the double four-brane bigravity models. In section four we generalize the quasi-localized 
and crystalline constructions in six dimensions. Finally, we summarize the six dimensional 
cases and conclude. 

2 Single brane models in six dimensions 

At first we will discuss the single brane solutions in six dimensions to get an insight on 
the more complicated multigravity configurations. In the following subsection we will re- 
view the minimal model where the bulk cosmological constant is homogeneous and then 
a more generalized model with arbitrary cosmological constant components along the four 
dimensional and the extra dimensional directions. In the following we will implicitly assume 
orbifolding around the four-brane. 

2.1 The minimal single brane model 

The simplest single brane model consists of a four-brane embedded in six dimensional AdS 
space [Q . One of the longitudinal dimensions of the four brane is compactified to a Planck 
length radius R while the dimensions transverse to the four-brane compact dimension are 



3 We would like to thank Panagiota Kanti and Luigi Pilo for bringing this issue to our attention. 
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Figure 1: The minimal single four-brane model warp factors a(p) and 7(p). 

infinite. The most general spherically symmetric ansatz that one can write down in six 
dimensions and which preserves four dimensional Poincare invariance is[]: 

ds 2 = a(p)rj fll/ dx fl dx u + dp 2 + ^{p)d9 2 (1) 



where 9 is the compactified dimension with range [0, 2tt] and p is the infinite radial dimen- 
sion. 

It is straightforward to solve the Einstein equations for the bulk energy momentum 



tensor T 



(B) N 
M 



-A5fj with the brane contribution: 



T 



(br) N 
M 



S(P) 



vs; 



(2) 



V 



From the form of the Einstein equations given in the Appendix we see that the solution 
for the two warp factors is: 



0(p) 



-kp 



7 (p) = R 2 e ~ k P 



(3) 



with k 2 = — 10A/4 , where M is the six dimensional fundamental scale and the arbitrary 
integration constant R is just the radius of the four-brane (see fig •[[]). The Einstein equations 
require the usual fine tuning between the bulk cosmological constant and the tension of the 
four brane: 



5k 



(4) 



6 Note that in this paper we use different metric signature and different definition of the fundamental 
scale from [1421 . 



Let us note at this point that in [42] this fine tuning was absent because a smooth local 
defect was considered instead of a four-brane. In this case, the fine tuning emerges between 
the different components of the defect energy momentum tensor. The physics of the four- 
brane idealization and the one of the defect model is the same. 

The four dimensional KK decomposition can be carried out as usual (see [S3| for a 
detailed fluctuation analysis) by considering the following graviton perturbations: 

ds 2 = a(p) [ripv + /v(P> #7 x )\ dx^dx u + dp 2 + 'j(p)d9 2 (5) 



Here, as well as throughout this paper, we have ignored the modulus associated to the 
radius R of the four brane. This will be a massless scalar and in order not to give rise to an 
additional long range force it should be given mass through some stabilization mechanism. 

We expand the graviton perturbations in a complete set of radial eigenfunctions and 
Fourier angular modes: 

V(P,M) = Y,^,i){p)e lW h^ l \x) (6) 

n.l 



The differential equation for the radial wavefunctions <fi is: 



V 



m 



R? 



? k p, 







(7) 



with normalization J dpa^^^n = 5 mn . We can convert this equation to a two dimen- 
sional Schrodinger-like equation by the following redefinitions: 

, * = a 3/ U (8) 



§W 



so that 



■2J=d,(yfid,*)+V eff * = T * 



V f 



21k 2 



eff\ 



32 (f + !)• 
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4-T'W (9) 



From the form of the potential (see fig.0) we can easily deduce that the angular exci- 
tations spectrum will consist of continuum starting from a gap of the order j^ and thus 
can be safely ignored. For the s-wave (/ = 0) there is a normalizable zero mode which is 
a constant in the p coordinate, i.e. <po = const. The KK tower for the s-waves will again 
form a continuum but this time gapless with wavefunctions given by |£| : 



N m e* kp 
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Figure 2: The form of the effective potential V e jf for different angular quantum numbers I 
for the minimal single four-brane model. The 5-function in the origin is omitted. 



The correction to the Newton's law due to the s-modes can be easily calculated and one 



finds that it is more suppressed than the five dimensional Randall-Sundrum case |42],[54|,[55 
In six dimensions the correction reads: 



AV 



0{M% 



ii: 



pi j 



2.2 The generalized single brane model 

We can now try to find more general single four-brane solutions by relaxing the requirement 
of homogeneity of the bulk and brane energy momentum tensors. In contrast to the five 
dimensional case, this is possible in six dimensions because there are two independent 
functions in the metric (see Appendix). We will consider the following inhomogeneous bulk 
energy-momentum tensor: 



T 



(B) N 

M 



Ao# 



A, 



(12) 



A, 



and allow for an inhomogeneous brane tension of the form: 



T 



(br) N 
M 



-Kp) 







(13) 



Va 



This inhomogeneity can be due to different contributions to the Casimir energy in the 
different directions 0,0 (see also p5| ) or due to a background three-form gauge field 
with non-zero field strength [HO]. 



2A fl 



If we define the parameter a = -^- £ , then the one four-brane solutions for the warp 
factors are: 



*(P) 



-kp 



j(p) = Re * 



|(5a-6)p 



(14) 



with k 2 = — tqttz , where again R is the radius of the four-brane (see fig.0). The flatness of 
the brane is achieved by the following fine tunings of the the parameters of the bulk and 
brane energy momentum tensors: 

8A« 



V ft 



5k 



_5« + 6 

VQ — T7, V0 

16 



5a 2 + 12 A 



(15) 



Let us note that the two components of the brane tension have the same sign as long as 



a > — |. For a > | the internal space is shrinking as in |4^], for a < I the internal space 
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Figure 3: The generalized single four-brane warp factors a(p) and 'j(p) for different values 
of the parameter a. 
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Figure 4: The form of the effective potential V e ff for a > 2 for different angular quantum 
numbers I. The 5- function in the origin is omitted. 



is growing and for a = | the internal space is a cylinder (see fig.0). 

Proceeding with the usual KK decomposition (see ||54|| ) we find the following equation 
for the radial wavefunction d>: 



J" - ^p%' + (m 2 - ^e^" 2 *' ) e k ^ = 



(16) 



with normalization j dpa y/^y (j) m (j) n = 5 mn . With the wavefunction redefinition and the 
coordinate change discussed in the previous subsection, we obtain the following effective 




Figure 5: The form of the effective potential V e ff for a < 2 for different angular quantum 
numbers I. The S- function in the origin is omitted. 
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potential for the two dimensional Schrodinger-like equation: 

T . , , 3(5« + 4)A: 2 I 2 fkz \& a -*> 3k .. . . 1P . 

For the s-wave there is a normalizable zero mode when a > — ■= which is a constant 

5 

in the p coordinate, i.e. 0o = const. The KK tower for the s-waves will again form a 
continuum and their wavefunctions are given by: 



N m eZ kp 



'"(tMM-^tWt' 1 ' 



where z/ = |(a + 2). The correction to the Newton's law due to the s- modes can be easily 
calculated and one finds that it is negligible as long asa> — | [54] : 



AV =-o(mr^ (19) 

As far as the angular excitations I ^ are concerned, in the case where a > 2 the 
potential diverges at infinity (see fig. |3j), so the spectrum will be discrete starting from a 

scale bigger than -jU- and thus can be safely ignored. On the other hand, the potential for 
the case a < 2 asymptotically vanishes (see fig. [5]) and so the spectrum will be gapless con- 
tinuum. One should examine the contribution of the angular excitations to the Newtonian 
potential in this case. 

3 Bigravity in six dimensions 

The flat single four-brane models considered in the previous sections have the characteristic 
that gravity is localized on the brane in the same way as in the five dimensional analogue 



55| . In this section we will show how we can construct realistic multi-localization scenarios 
for gravity by sewing two single four-brane solutions together. We will present two explicit 
examples of compact double four-brane models where bigravity is realized. Firstly, we will 
consider a double four-brane bigravity model which in general contains an additional three- 
brane associated with the existence of a conical singularity. This can be done by allowing 
for an inhomogeneous four-brane tension. In the following subsection we demonstrate how 
one can avoid the conical singularity by also allowing for inhomogeneous bulk cosmological 

10 



constant. In the following we will consider the orbifold geometry S 1 jZi in which the two 
four-branes lie on the fixed points. 

3.1 The conifold model 

We are interested in a double four-brane generalization of the minimal single four-brane 
model. In order to achieve this, we still assume a homogeneous bulk energy momentum 

tensor T M = — A5j^ but we do not impose any constraints on the four-brane tension: 

T {b V -S{p) [ (20) 

V e 

The Einstein equations for a(p) and 7(p) in this case give the following solutions for the 
warp factors: 

cosh 4 / 5 [ jjfc(p - p )] cosh 6 / 5 [jfcpp] sinh 2 [f fc(p - Po)] 

[P) ~ cosh 4 / 5 [ffcpo] ' 7lPJ sinh 2 [ffcp ] cosh 6 / 5 [f fc(p - p )] l j 

with k 2 = — 10M 4 , where we have normalized cr(0) = 1 and 7(0) = R 2 (see fig. |^). From 
the above relations it is obvious that both <r(p) and 'y(p) have a bounce form. However, 
we note that 7(po) = 0, that is, 'j(p) is vanishing at the minimum of the warp factor cr(p). 
This is a general characteristic of the solutions even when the branes are non-flat. From 



eq. (|A.3|) and eq.( [A.4|) we can easily show that 7(p) = C Jfl (where C is an integration 



constant) which implies that whenever we have a bounce in the warp factor the function 
7(p) will develop a zero. 

In order to examine the nature of this singularity we study the form of the metric at 
the vicinity of the point p = po along the lines of [p3,|B0]. Taking in account that in this 



limit we have a(p) — > b 2 and 7(p) — > /3 2 (p — po) 2 the metric becomes: 

ds 2 = b 2 r ]llu dx> l dx' / + dp 2 + /3 2 (p-p ) 2 d6 2 (22) 

where b 2 = cosh -4 / 5 [ffcpo] 

and 1 = 



_ 25fc 2 J R 2 cosh 6 / 5 (|fcp ) 
16 sinh 2 (|fcp ) 



From the form of the metric, it is clear that for general values of the j3 parameter there 
will be a conical singularity with a corresponding deficit angle 5 = 2ir(l — (3). The existence 

11 





Figure 6: The conifold bigravity model warp factors a(p) and 'y(p) for the symmetric 
configuration where L = 2p . Both the a(p) and 7(p) have a bounce at po- However, 7(p) 
vanishes at this point. This point in general corresponds to a conical singularity. 



of this conifold singularity is connected to the presence of a 3-brane at p = po- In order 
to find its tension one has to carefully examine the Einstein tensor at the vicinity of the 
singularity. For this reason we write the metric for the internal manifold in the conformally 
flat form: 



ds 2 = tfrj^dx^dx" + f{r){dr 2 + r 2 d9 2 ) 



(23) 



with /(r) = r 2 ^ 1 - ) and p — po = P l ^ ■ In these coordinates it is easy to see how the three 
brane appears at the conifold point. The Einstein tensor can be calculated for p — > p : 



R 



MN 



-GmnR 



V 2 log(/(r)) h 2 _ 
2 fir) U 'l^ v 



(24) 



where V 2 is the flat two dimensional Laplacian. Now, given that V 2 log(r) = 2-k5^ 2 \y) and 



12 



by comparing with: 



Run ~ \GmnR = -^^G^ M 5 v N 5{r) (25) 



where 5^ 2 \r) = -$- (see Appendix of @), we find that the tension of the 3-brane is: 

V 3 = 4(1 - /?)M 4 = ^—5 (26) 

Thus, if there is angle deficit 5 > ((3 < 1) the tension of the brane is positive, whereas 
if there is angle excess 5 < (j3 > 1) the tension of the brane is negative. At the critical 
value (3 = there is no conical singularity at all and we have a situation where two locally 
flat spaces touch each other at one point. 

In the following we will concentrate on the symmetric case, that is, the four-branes will 
be considered placed at symmetric points with respect to the minimum of the warp factor 
p , ie. L — 2pQ. For the previous solution to be consistent the tensions of the four-branes 
for the symmetric configuration must satisfy: 

3„ 8A 2 1 



Vg = — — — tanh 
ok 



2 kp0 



V °=8 V ' + WV. (27) 



The brane tensions for both branes are identical. Thus, the above construction consists of 
two positive tension four-branes placed at the end of the compact space and an intermediate 
three-brane due to the conifold singularity with tension depending on the parameters of the 
model. In the limit po —* °° we correctly obtain two identical minimal single four-brane 
models for the case where 5 — > 2-k {(3 —* 0). 

The differential equation for the radial wavefunction <fi of the graviton excitations is: 



with normalization J Q Po dpo ^<\) m <\) n = 5. 



•mn ■ 



There is an obvious normalizable zero mode with 0o = const, and a tower of discrete 
KK states. It is easy to demonstrate that for 5 > ((3 < 1) this model has an ultralight 



state, leading to a theory of bigravity without negative tension branes. This follows [|61 
considering the limit po ~~ * oo when the theory describes two identical single four-brane 
models each of which has a massless graviton. For finite po only their symmetric combination 

13 
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Figure 7: The form of the effective potential V e ff for different angular quantum numbers 
/ for the conifold bigravity model in the symmetric configuration. The ^-functions in the 
positions of the four-branes are omitted. 



remains massless but their antisymmetric combination acquires an anomalously small mass, 
because the difference in the modulus of the wavefunctions between the massive and the 
massless states is significant only near po where the wavefunctions are exponentially smallQ 

We can transform this differential equation to a two dimensional Schrodinger-like one 
by the coordinate change ^ = <t _1//2 = g(z) and the usual wavefunction redefinition. Then 
the effective potential reads: 



V, 



eff 



15 fd z g 
9 



3d 2 z g 3 d z gd zl | I 2 
4 g 8 37 2 ig 2 



(29) 



We cannot write an explicit analytic formula for of the above potential in the z- 

coordinates because the coordinate transformation is not invertible analytically. However, 

since the transformation is monotonic, we can easily sketch the form of the potential (see 

fig.[7|) by calculating it in the p coordinates. From this procedure we see that the potential 

for the s-wave is finite at the conical singularity po but has a divergence for all the angular 

excitations with / ^ 0. This means that only the s-wave excitations with communicate the 

There is, however, a mismatch of the degrees of freedom of the massive and massless graviton, which 
may not be phcnomcnologically dangerous, but still raises the issue of locality in the infinite four-brane 
separation limit. 
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two parts of the conifold and the other excitations will be confined in the two semicones. 
The spectrum of KK states in this case will be discrete. 

The construction may be readily modified to obtain a six-dimensional analogue of the 



locally localized model ||19|| . This is done by considering the asymmetric situation in which 
one of the fixed point four-branes is moved to infinity (in the asymmetric case the tensions of 
the two four-branes are different). Then gravity on the four-brane at p = will be mediated 
by only the ultralight state since the graviton zero mode will not be normalizable. 

3.2 The non-singular model 

The appearance of a conifold singularity can be avoided if instead of using the minimal 
single four-brane solution to build the double four-brane solution, one uses the generalized 
single four-brane model. By allowing for an inhomogeneous bulk cosmological constant one 
can arrange that the function 7(p) does not develop a zero. In this case the (6, 8) and (p, p) 
part of Einstein equations in the bulk, ( |A.3| ) and (|A.4| ), have solutions: 

cosh 4 / 5 [jkjp - p )] sinh" [\k{p - go)] 

a (P> = u4/5r5; 1 ' ™ = A ,6/5 5,, 77 30 ) 

cosh 4/5 [ffcpo] cosh 6/5 [\k{p - po)] 

where A is a constant and a = -z- 2 -. The (p, v) component of the Einstein equations (|A.5|) , 

however, restricts the possible values of a to a = 2 or a = 0. The case of a = 2 corresponds 
to the previous conifold model. For the case a = 0, which corresponds to the choice A p = 0, 
the solution is: 

cosh 4 / 5 [jkjp - po)] n2 cosh 6 / 5 [f kp ] 

a ( p ) = ,4/5r5> i ' 7 W = R — ,6/5 rsT7 ^T 31 

cosh 4/5 [ffc/oo] cosh b/5 [\k{p - po)] 

with k 2 = — YoM*' where we have normalized cr(0) = 1 and 7(0) = R 2 (see fig. ||). From 
the above the absence of any singularity is obvious since ~f(p) does not vanish at any finite 
value of p. 

The four branes will appear as ^-function singularities at the points where we will cut 
the space in the p direction. In the following we will concentrate on the symmetric case, 
that is, the branes will be considered placed in symmetric points in respect to the extremum 
of the warp factor p , i.e. L = 2p . In order for the above configuration to be realized, the 
components of the four-brane tensions and the bulk cosmological constant must be tuned 
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Figure 8: The non-singular bigravity model warp factors a(p) and j(p) for the symmetric 
configuration where L = 2p . The <r(p) warp factor has exactly the same form as in the 
conifold model. The 'y(p) warp factor however now has an inverse bounce form. In this 
case 7(p) does not vanish anywhere and thus the model is free of singularities. 



to give: 



8A , 
Vg = tanh 

5k 



-kpo 



In 



;Vg 



Ar 



:A fi 



(32) 



The brane tensions for both branes are identical. Thus, we have constructed a compact 
model with two positive tension flat branes with a bounce in the warp factors. Note the 
extra condition relating the different components of the bulk cosmological constant; this 
inhomogeneity in the bulk cosmological constant is what allows us to have a non-singular 
solution. In the limit pn ~ > °° we obtain two identical generalized single four-brane a = 
models. 
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v< 



eff 




2f) 



2z 



Figure 9: The form of the effective potential V e ff for different angular quantum numbers 
/ for the non-singular bigravity model in the symmetric configuration. The (^-functions in 
the positions of the four-branes are omitted. 



The differential equation for the radial wavefunction <fi of the graviton excitations is 
given by (|28| ) for the relevant a(p) and 'j(p) functions. Although it is difficult to find an 
analytical solution of the corresponding differential equation for the KK states, the existence 
of a light special KK is assured from the locality arguments discussed above ||6l |. Since 
at the infinite brane separation limit we recover two identical generalized single four-brane 
models (for a = 0) where each of them supports a massless zero mode, at a finite separation 
configuration the one of the zero modes will become the special light KK state. 

By making an appropriate change of coordinates and a redefinition of the wavefunction 
we can bring the previous differential equation in a two dimensional Schrodinger form with 
the effective potential ( |29"D that is plotted in fig|| Note that in this case the angular 
excitations feel the whole region between the two positive branes and thus communicate 
physics between them. The spectrum of KK states in this case will again be discrete. 



4 Multigravity in six dimensions 

In this section we will show how we can obtain the analogues of the GRS model of quasi- 
localized gravity and of the crystal universe in six dimensions, in both cases without the 



17 




Figure 10: The a = quasi-localized model warp factors cr(p) and 7(p). The bulk region 
for p > p b is flat. 



need of introducing moving negative tension branes. 



4.1 Quasi-localized gravity 



As in the bigravity models we have two possible ways of realizing a quasi-localized scenario 
in six dimensions. The first model can be built from the a = bigravity model, by cutting 
the space at a point p b on the right of the position of the bounce and sewing it to flat space 



on the right (see fig.|10]). In that case we will have a system of two positive tension branes 
and by appropriately tuning the different components of the brane tension of the second 
brane we can find solutions where the warp factor is constant in the flat bulk region. The 
fine tunings which achieve this are: 



V t 



(2) 



4A , 

tanh 

5k 



Mpb - po) 



V r 



(2) 



IvP 



(33) 
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Figure 11: The a = 2 quasi- localized model warp factors a(p) and 7(p). The bulk region 
for p > pb is flat. 

From the above formulas it is obvious that if we paste the flat space at the position of 
the bounce, i.e pt = Po then the second brane is tensionless. If we relax these fine tunings 
we can find solutions where the warp factors in the flat bulk region are non trivial, but we 
will not discuss this possibility here. 

The above construction will obviously have no normalizable zero mode due to the infinite 
volume of the system and the KK spectrum will be continuous. However, as in the five 
dimensional case we expect that the low part of the KK tower will have a "resonance" - 
like coupling to matter on the "+" brane at p — 0. This will mediate normal gravity at 
intermediate distances but will change the nature of the gravitational law at ultralarge 
scales. 

The second model of quasi-localized gravity can be built in a similar fashion from the 
a = 2 bigravity model (see fig.jll]). Again we obtain a system of two positive tension branes 



which will have a constant warp factor in the flat bulk region when the following fine tunings 
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Figure 12: The a = crystalline model warp factors <x(p) and 7(p). 



are demanded: 



v; 



(2) 



4A 
~5k 



tanh 



Mpb ~ Po) 



V f 



(2) 



8 e 



8A 2 



5/c 2 y 



(2) 



(34) 



From the above formulas it is obvious that if we paste the flat space at the position of the 

(2) (2) 

bounce, i.e pb = po then the tension V@ vanishes whereas the tension V^ diverges. This 
limit is not well behaved since the internal space becomes singular in a whole line. At this 
limit obviously classical gravity breaks down and one would expect that quantum gravity 
corrections would resolve the singularity. However, one should check if these quantum 
gravity corrections affect the solution in the vicinity of the four-brane at p = 0, e.g. by 
examining the R 2 correction to the Einstein actionQ 

The above construction will again have no normalizable zero mode and the KK spectrum 
will be continuous. A "resonance" -like coupling of the KK states to matter on the "+" brane 
8 We would like to thank Tony Gherghetta for this comment. 
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Figure 13: The a = 2 crystalline model warp factors a(p) and 7(p). 

at p = is again expected which will give normal four dimensional gravity at intermediate 
distances but will change the nature of the gravitational law at ultralarge scales. 



4.2 The crystal universe model 

Another obvious generalization of the five dimensional multigravity models is the one of 
the crystalline brane model. All that one has to do is to paste an infinite array of bigravity 
models. We again take this constructions in two copies. One for a = (fig.|12|) and one 
for a = 2 (fig.|T3"D. A band structure is again expected as in the five dimensional case and 
the width of the first band will be exponentially smaller than the one of the first forbidden 
zone. As far as the phenomenology is concerned, these models will generate normal gravity 
at intermediate scales as the first band will behave as an effective zero mode, whereas 
at ultralarge scales will have the same change of the gravitational law as the one of the 
corresponding quasi-localized models. 
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5 Discussion and conclusions 

In this paper we have constructed, for the first time, flat brane theories which can lead to 
multigravity models and their associated modifications of gravity at large distances without 
introducing moving negative tension branes. The constructions are made possible by going 
to a six dimensional theory. In five dimensions with flat branes, the presence of a bounce 
of the warp factor was linked to the violation of the weaker energy condition. This is not 
true, however, in six dimensions as we will show in this section. In five dimensions with 
the metric: 

ds 2 = e - A{p) Vl , u dx»dx u + dp 2 (35) 

one can readily show that the weaker energy condition requires that: 

A" > (36) 

which is violated at the position of the moving negative tension branes. In the case that 
the branes were AdS one could have a bounce (without having negative tension branes) 
and still satisfy the weaker energy condition because the above relation is modified to [TJJ : 

A" > -2H 2 e A (37) 

However, such models do not lead to modifications of gravity at large distances and moreover 
the remnant negative cosmological constant is in conflict with current observations. 

In the six dimensional case with metric: 

ds 2 = e- A ^ri^dx^dx u + dp 2 + e~ B{p) d6 2 (38) 

from the weaker energy condition one finds two inequalities which can be cast into the 
following relation: 

h B ') 2 - -B" - -A'B' < A" < ~-(B') 2 + B" - -A'B' + 2{A'f (39) 

6 3 6 2 2 

These inequalities hold for the a = 2 models as long as the three-branes sitting on 
the conical singularities have positive tension. It is violated, however, in the a = models 
everywhere in the bulk. This shows that in six dimensional flat brane models, the presence of 
bounces of the warp factor is not necessarily linked to whether the weaker energy condition 
is satisfied or not. It is not yet clear if the violation of the weaker energy condition in the 
a = models is a sign of a possible instability. 
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There is an interesting feature of the conifold model that emerges if we consider Eu- 
clidean four dimensional space and Wick rotate the 9 coordinate (8 —>■ it) decompactifying 
it at the same time. In the case where the deficit angle is zero, the metric in the vicinity of 
p = Po will look like: 

ds 2 = b 2 5^dx^dx v + dp 2 - (p - p ) 2 dt 2 (40) 

This is Rindler space which is deformed away of p = po an d has sources (the four-branes) at 
finite proper distances from p = p . The picture that we obtain in Kruskal coordinates for 
this case is reminiscent of the one of a Schwarzschild black hole if one rotates the Penrose 
diagram at a right angle. 

In summary, in this paper we discussed how we can obtain viable multigravity models 
in six dimensions. The models we considered involved flat positive tension branes (except 
the case when the conifold models had angle excess). We showed how we can obtain the 
double brane bigravity model, the quasi-localized model and the crystalline model, each of 
which in two types. One of them is singularity free while the other has conical singularities 
in the bulk corresponding to three-branes. 

Several issues are still open regarding these constructions. Firstly, a careful treatment 
of the angular excitations in the case of a = should be carried out to examine the 
nature of gravity on the four-branes. Furthermore, one should examine the moduli of the 



system, the radions and the dilaton, and calculate their mass (see [[T3|,[62[] for the warped 
five dimensional case and |H[] for the warped six dimensional case). It is quite probable 
that the modulus corresponding to po is massive as in the five dimensional AdS± branes case 
and so the system is self-stabilized. Moreover, a discussion of multigravity in a cosmological 
setting is yet to be developed (for a discussion about cosmology in brane world models see 



for example [63—70]). Finally, the non-trivial way that locality is preserved by these model 
should be studied. 
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Appendix 

In this Appendix we provide the Einstein equations for the class of models that we 
considered where the metric is: 

ds 2 = a(p)rj fll/ dx f "dx" + dp 2 + ^{p)d9 2 (A.l) 

Ignoring the three-branes located on the conical singularities, the Einstein equations are 
generally written as: 

G MN = -^i( T MN + T MN> ( A - 2 ) 

where T^h is the bulk energy momentum tensor of the general form ([12]) and T^iv the 
four-brane energy momentum tensor of the form fll3D . 

In the absence of a four dimensional cosmological constant, the (9, 9) component of the 
above equation is: 

a" 1 fo'V A Vi 



2 T + 5 7 =-«F"^("-"- ) (A ' 3) 



The (p, p) component is: 

3 /rr'\ 2 n-'V A 

(A.4) 



3 fa'Y a'i _ A p 



2\a J (77 4M 4 

Finally, the (p, v) component is: 

3a" 3a'i 1/yV ,1-f Ao Vj 

These equations may be compared with the ones of the five dimensional case where the 
metric: 

ds 2 = a(p)rj^dx^dx u + dp 2 (A.6) 



(A.7) 



gives rise to the (p, p) component: 








l®'- 


A 
4M 3 


and the (p, v) component: 






3 a" 


A 


V i 



2 a 4M3 AM* 8{p - pi) (A ' 8) 

The extra freedom that we have in the six dimensional case is apparent. 
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